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^ . FUSE spectra of star-forming regions in nearby galaxies are compared to 

^ . composite spectra of Lyman-break galaxies (LBGs), binned by strength of Ly- 

man a emission and by mid-UV luminosity. Several far-UV spectral features, 
^ . including lines dominated by stellar wind and by photosplieric components, are 

^ . very sensitive to stellar abundances. Their measurement in Lyman-break galax- 

\Q , ies is compromised by the strong interstellar absorption features, allowing in 

some cases only upper limits. The derived C and N abundances in the LBGs 
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\0 . are no higher than half solar (scaled to oxygen abundance for comparison with 

^ , emission-line analyses), independent of the strength of Lyman a emission. P V 

'"q_^! absorption indicates abundances as low as 0.1 solar, with an upper limit near 

0.4 solar in the reddest and weakest-emission galaxies. Unresolved interstellar 
absorption components would further lower the derived abundances. Trends of 
d '. line strength, and derived abundances, are stronger with mid-UV luminosity than 

with Lyman-a strength. H2 absorption in the Lyman and Werner bands is very 
^ , weak in the LBGs. Template H2 absorption spectra convolved to appropriate 

^ ! resolution show that strict upper limits N(H2) < 10^^ cm~^ apply in all cases, 

with more stringent values appropriate for the stronger-emission composites and 
for mixes of H2 level populations like those on Milky Way sight lines. Since the 
UV-bright regions are likely to be widespread in these galaxies, these results rule 
out massive diffuse reservoirs of H2, and suggest that the dust /gas ratio is already 
fairly large at z ~ 3. 

Subject headings: galaxies: evolution - galaxies: abundances — ultraviolet: 
galaxies 



^Based on observations made with the NASA-CNES-CSA Far Uhraviolet Spectroscopic Explorer. FUSE 
is operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985. 
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1. Introduction 

Massive stars play crucial roles in galaxy evolution, quite aside from dominating many 
of the observed properties of star-forming systems. Their mass loss and energy input are 
major factors in the history of the interstellar medium. The far-ultraviolet range, between 
the Lyman limit and Lyman a, shows us only these stars, giving us a direct window to the 
ongoing star formation wherever the extinction is low enough. The far-UV range provides 
an embarrassment of riches for tracing stellar winds and the interstellar medium in both 
atomic and molecular form, inchiding atomic species and ionization stages not represented 
at longer wavelengths. Both extinction and the history of star formation make only subtle 
differences in the observed far-UV spectrum, in each case because the parameter range over 
which stars contribute to the far-UV light is rather limited. 

The rich information content of FUSE spectra has revolutionized our understanding 
of the UV universe, particularly including star-forming systems. In this paper, I compare 
properties of local star-forming galaxies, whose chemical compositions are well established 
from emission-line spectroscopy, to those of Lyman-break galaxies aX z ^ 3. This compar- 
ison allows estimates of the metallicity in the massive stars, and sets interesting limits on 
absorption from surrounding molecular gas. 

2. Data 

2.1. FUSE spectra 

The local reference data are FUSE observations of luminous star-forming regions in 
nearby galaxies, described by Keel, Holberg, & Treuthardt (2004). These consist of the 
four UV-brightest giant H II regions in M33 (NGC 588, 592, 595, and 604), a bright outer 
H II region in MlOl (NGC 5461), the starburst nucleus of NGC 7714, and the summed 
spectrum of I Zw 18 produced by Aloisi et al. (2003). These spectra were all obtained 
with the large 30" aperture, encompassing essentially the entire young populations with 
the goal of comparison with integrated spectra of distant galaxies. Exposure times for the 
new observations were 4000-7000 seconds, while I Zw 18 was observed for a total of 95,000 
seconds. Since NGC 604 is an order of magnitude brighter than the other objects in the 
FUV band, its spectral signal-to-noise ratio is the best of the whole sample, and the data 
can be evaluated at higher resolution (0.1 A ). As detailed by Keel et al.. Galactic H2 hnes 
were removed by fitting families of Gaussian or Voigt profiles (for the stronger transitions) to 
the various H2 multiplets, with line identifications from the H2OOLS package by McCandliss 
(2003). Several atomic ISM lines from the Milky Way and high- velocity clouds were similarly 
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excised. 

Fig. 1 identifies important interstellar and wind features in the far-UV spectral range 
on the FUSE spectra of NGC 604 and I Zw 18. These objects play special roles in our 
understanding of hot stellar populations in the far-UV (and indeed at other wavelengths 
as well). NGC 604 is by far the most luminous star- forming region in the nearby spiral 
NGC 598 (M33), comparably powerful to the 30 Doradus complex and perhaps NGC 3603 
in our own Galaxy. However, it does not contain a central super star cluster, but is instead 
a collection of smaller and less dense subsystems, whose star-forming history featured a 
strong burst about 3 Myr ago (Gonzalez Delgado & Perez 2000)). Active star formation 
has proceeded long enough in NGC 604 to generate an overall bubble geometry in the 
surrounding H II region, one so extensive that it probably contributes to the low limits on 
foreground H2 through association or displacement of the surrounding disk ISM (Bosch et 
al. 2002). The optical emission lines yield an oxygen abundance O/H^O.66 solar (Vilchez et 
al. 1988). In contrast to this relatively metal-rich object, I Zw 18 remains the galaxy with 
the lowest known emission-line abundances (0/H only 0.05 solar), more than 30 years after 
the measurements of Searle & Sargent (1972) showed its exceptional nature. Studies of its 
colors and resolved stars have often centered on asking whether there is evidence for any star 
formation before the recent burst which dominates its luminosity. The question is whether I 
Zw 18 has only recently come into view as a galaxy of stars, or whether this is simply a very 
extreme example of a starburst in a previously faint but gas-rich dwarf system. Whatever 
its previous history, at such low metallicities, the abundance pattern is dominated by only a 
sparse history of enrichment, making it potentially a useful analog to chronologically young 
galaxies seen at high redshifts. 

Table 1 lists some of the most prominent ionic absorption species seen in the far-UV, 
identified as to main source (photospheres, winds, or interstellar absorption). Many have 
multiple sources (particularly from either winds or the ISM) . The wavelength pattern shows 
the blending of several otherwise promising sets of features which means that only about 
five lines are both strong enough and isolated enough to examine in the (low-dispersion) 
data available for high-redshift galaxies. The NGC 604 data in particular suggest, via the 
wavelength centroid, that P V at 1128 Ais significantly stronger than the adjacent Si IV 
feature (which is why the former identification is used here). 

Although some velocity confusion exists between M33 absorption and a foreground high- 
velocity clouds, some interstellar features seen in the NGC 604 spectrum appear to arise in 
the surrounding disk gas. Likewise, there are some interstellar lines from the ISM in I Zw 
18 (although most of the prominent ones seen in Fig. 1 arise on our own galaxy). These 
spectra provide only upper limits to intrinsic H2 absorption, at column densities less than 
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10^* — 10^^ cm~^, in agreement with the analysis of the same M33 spectra by Bluhm et al. 
(2003). The roles of both kinds of ISM absorption are important in comparison with LBG 
spectra, where both aperture and column-density effects enter. 

The excellent spectral resolution of the FUSE data helps clarify some potential issues of 
line blending for actual hot stellar populations. For example, they show that the P V A1128 
feature, which for most stellar types is photosphcric. is much stronger than the adjacent Si 
IV 1122, 1128 pair or Si III A1113. This contrasts with the situation in the LBGs, where 
interstellar absorption lines are much stronger and in some cases overlap with diagnostic 
abundance-sensitive lines. These interstellar features, generally showing strong outflows, are 
the major difference in comparing the local star-forming objects with Lyman-break galaxies. 
The column density of foreground gas is an extensive property, which naturally tends to 
increase with galaxy scale, making this one respect in which individual luminous star-forming 
regions may not be adequate analogs of Lyman-break galaxies. 



2.2. Lyman-break galaxy spectra 

The FUSE spectra are compared to the group-composite spectra of Lyman-break galax- 
ies (LBGs) discussed by Shapley et al. (2003), spanning the emitted wavelength range 
912-1200 A . The major variation among spectra of LBGs is in the net strength of Lyman 
a emission, which has variously been attributed to effects of metallicity, geometry, and the 
importance of global winds. As shown in Table 3 of Shapley et al., the sequence of Lyman 
a emission is also statistically a ranking by continuum color (and hence reddening). The 
composite spectra include 811 galaxies at 2; 3, summed in quartiles (roughly 200 galaxies 
each) of Lyman a equivalent width. The sum not only increases the signal-to-noisc ratio of 
spectral features in common, but is crucial in averaging out the effects of absorption in the 
Lyman a forest. The quahty of these spectra varies, since there is a rough inverse relation 
between Lyman a EW and continuum luminosity. 

I also consider a set of composite spectra grouped into three bins by mid-UV luminosity. 
Given the present-epoch increase of metallicity with luminosity, this subsampling might trace 
chemical history more purely than the division on Lyman a strength, which may be driven 
by such transient properties as global winds or dust distribution. However, this division 
is subject to dust effects, so that one may fairly wonder whether the UV luminosity scales 
strongly enough with star-formation rate, less alone stellar mass, to reflect any systematic 
chemical behavior. Given that the UV luminosity must be high enough for selection as 
an LBG, we might at least test for luminosity-linked effects. Nearby galaxies studied with 
lUE do show evidence that the mid-UV luminosity has at least an envelope, and perhaps 
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a broad correlation, with emission- line abundances (Calzetti et al. 1994, Keel et al. 2002). 
These luminosity composites are derived from the subset of LBGs near the peak of the color- 
selection function with redshift, 2.9 < z < 3.1, and will unavoidably reflect spectroscopic 
biases (most notably the easier measurement of redshifts for emission-line objects at the 
faintest levels, discussed by Shapley et al. 2003). Absolute magnitudes are evaluated in 
the ABj, system near 1700 A , using a cosmology with Hq = 71 km s~^ Mpc~^, Vim = 0.3, 
= 0.7. The faintest bin, spanning from Muy=-20.55 to -19.52, includes 90 galaxies. The 
intermediate-luminosity bin samples from M[/y=-21.02 to -20.56 with 93 galaxies, and the 
brightest group includes 92 galaxies from M[/y=-22.58 to -21.03. Average Muv values are 
-20.25, -20.79, and -21.44 for these subgroups. As shown in Fig. 3, some strong absorption 
features do exhibit systematic trends with UV luminosity. The O VI doublet and both 
N II and N III (the latter possibly blended with interstellar O I) all strengthen for more 
UV-luminous galaxies, results which suggest that there may indeed be composition changes 
along this sequence. The signal-to-noise ratios of the composite spectra are naturally rather 
different; there are only rather loose upper limits for the bluest lines (AA955, 991) for the 
fainter groups. 

The effective resolution of these spectra is limited both by the uncertain placement 
of the galaxies within fairly wide slits and by uncertainty in the individual redshifts, since 
global winds and radiative-transfer effects can make redshifts measured from Lyman a or 
ISM lines alone differ from the actual systemic value. I estimate the resolution of each 
spectrum through Gaussian fits to absorption features dominated by interstellar components, 
particularly Fe II and Si II. Most such lines are in the mid-UV from 1260-1670 A ; several 
have at least weak P Cygni wind-emission components. The overall average FWHM of 
these lines is 3.3 A , with a trend for narrower profiles in the stronger-emission composites. 
However, some of these lines are weaker in the objects with strong Lyman a emission, so these 
measurements arc correspondingly less certain. The effective line width may be dominated in 
some galaxies by the width of absorption lines from outflowing gas, which reaches a detected 
span of 1300 km s"^ in at least MS1512-cB58 (Pettini et al. 2002). Such widths would apply 
to interstellar absorption but not to stellar features. 



3. StellEir Abundances 

In contrast to the nebular emission lines which are prominent in the optical regime, 
the continuum from massive stars and both emission and absorption lines from their winds 
dominate the UV, and particularly far-UV, spectral region of star-forming galaxies. Reliance 
on absorption features makes abundance analyses observationally challenging, since emission 
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lines may place much looser demands on the signal-to-noise ratio in the continuum. However, 
the results are astrophysically easier to interpret, since wc arc guaranteed to measure the 
composition of the massive stars themselves at the epoch observed, rather than some mixture 
of pristine and enriched interstellar material. The strong stcllar-wind lines have particular 
promise as metallicity indicators, since abundances enter both into the fraction of the wind 
column density from each element and into the overall mass-loss rate in radiatively-driven 
winds. This notion has been applied by Pettini et al. (2000) and Leitherer et al. (2001), the 
latter using a synthetic library of model spectra. Mehlert et al. (2002) used the C IV A1550 
doublet to derive rough abundances for a sample of LBGs. They found carbon abundances 
evolving in the mean from 0.16 to 0.42 solar over the range z — 3.2 — 2.3, but use of the 
C IV lines alone, at low resolution, is likely to be compromised by the role of absorption 
from outflowing material, and be affected by the star-forming history and geometry of the 
galaxies, through possible variations in ionization balance and relative amounts of emission 
and absorption. Use of the far-UV range can expand this approach to a wider range of 
elements and ionization states. The grounds for comparison are shown in Figs. 2 and 
3, which show the locations of wind, photospheric, and interstellar lines on the two sets of 
composite LEG spectra, in comparison to a similarly smoothed version of the NGC 604 FUSE 
spectrum. In particular, the LEG composites show stronger interstellar atomic absorption 
than any of the local FUSE targets, which for some features limits the accuracy of derived 
absorption equivalent widths for the stellar features. 

Keel et al. (2004) used the FUSE results to propose an empirical relation between 
the equivalent widths of far-UV photospheric and wind lines and the overall metallicity, 
parameterized through [0/H] as measured from the optical emission lines. These equivalent 
widths are with respect to local pseudocontinuum bands selected from the high-resolution 
FUSE spectra to be free of discrete absorption features, bands which were matched as nearly 
as possible in measuring the LEG composite spectra. Not all strong stcllar-wind lines in the 
far-UV are amenable to such analyses; for example, the components of the O VI doublet 
are usually thoroughly blended, with each other and with interstellar absorption of C II 
and Lyman /3. Useful relations exist for C II, C III, C IV, N II, N III, and P V. Among 
these, P V has a substantial photospheric contribution for many stellar types and has the 
most nearly linear EW-(0/H) behavior. As a group, all these relations should yield the 
same results only for galaxies with a long history of stellar nucleosynthesis and recycling, 
since various elements should enter the ISM and be incorporated into stars at quite different 
rates. In fact, one can use changes in the abundances derived from various species to seek 
signs of this history in high-redshift galaxies. Carbon and phosphorus are primary elements 
(although the yield may depend on stellar abundances, from the work of Maeder 1992), while 
the situation for nitrogen remains ambiguous at low metallicity. 
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Within the LBG composite spectra, the major changes with Lyman a strength are in 
the P V A1122 and N II A1084 lines, both of which weaken with stronger Ly a. Tables 2 and 3 
summarize the line properties of the LBG spectra, grouped by Ly a strength and by mid-UV 
luminosity. At this spectral resolution, absorption and emission components are seldom well 
separated for wind features, so net equivalent widths are listed. These should be robust to the 
resolution of the data as long as there are not issues of other spectral features contaminating 
the "continuum" windows. As in deriving the empirical EW-metaUicity relations, I have 
been guided by the FUSE spectra in selecting lines for which such blending issues are least 
important. 

As an extreme example of potential absorption blends, the Si II A 1526 line is nearly as 
strong as C IV wind absorption in most of the LBG composite spectra. Using large- aperture 
data for objects in the local sample from the lUE archive (as collected by Keel et al. 2004) to 
sample apertures comparable to the FUSE data, the ratio of equivalent widths between Si II 
and C IV absorption components falls in the range 0.38-0.50. In one LBG, MS 1512-cB58 at 
z = 2.728, gravitational Icnsing has allowed unusually detailed spectra to be obtained using 
Keck (Pettini et al. 2002) and the VLT (Savaglio, Panagia, & Padovani 2002), in an excellent 
example of Zwicky's (1937) "gravitational telescope". Far-UV lines in cB58 can suffer from 
confusion with the absorption in the Lyman a forest at this redshift, although upper limits 
can be set for emitted wavelengths longward of 1073 A for those lines which do not overlap 
an obvious lower-redshift Lyman a feature. However, these data are very informative about 
absorption hues longward of Lyman a. The Si II/G IV ratio of 0.48 (including the diluting 
effects of G IV emission, to mimic lower-resolution data) hes at the high end of what we see 
in local objects. 

The N II A1083 region in cB58 is more complex. Fairly broad absorption is present 
roughly centered on the expected wavelength of the (blended) N II triplet, with a strong 
narrow feature centered at Aq = 1080.0 A. This might be Fe II A1081.87 blueshifted by 

500 km s~^; if not, it would be the strongest observed feature in the Lyman a forest 
and this Fe II would be undetectably weak. This narrow, potentially interstellar component 
contributes a fraction 0.25 of the equivalent width of the whole complex, and would be 
thoroughly blended with the N II feature in our composite spectra. 

The equivalent widths in Tables 2 and 3 have been combined with the empirical relations 
from Keel et al. (2004) to estimate abundances for each of these species. For comparison with 
optical emission-line results, these relations use [0/H] as the independent variable. While 
convenient, it does leave the meaning of these abundances slightly opaque (being the abun- 
dance of that species typically encountered at the given oxygen abundance in present-day 
galaxies). Gomparing the multiple ions available for G and N is useful, since the ionization 
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state of winds will generally change with stellar abundances; the mean is listed for each as 
well as individual values. 

Effects of unresolved interstellar absorption hues will all increase the derived abun- 
dances, but the amplitude of the overestimate depends strongly on the mean velocity offset 
and width of this absorption component. A narrow component near the maximum stellar- 
wind absorption may have little effect since there is little light to absorb at that wavelength, 
while broad and blueshifted outflow proflles can demonstrably contribute a substantial part 
of the overall equivalent width. For resonance lines, this issue means that these derived abun- 
dances must strictly be considered as upper limits. This has been reinforced by the analysis 
of Crowther et al. (2006), in showing the importance of accounting for interstellar absorp- 
tion in the C IV transitions before deriving abundances to avoid a systematic overestimate. 
Once again, the scale of the system and its ISM can make a significant difference between 
spectroscopic properties of individual H II regions and entire galaxies (as seen for the LBG 
samples). In this comparison, the data for NGC 7714 and I Zw 18 encompass almost all their 
UV flux, while NGC 604 contributes nearly 20% of the entire far-UV flux from M33 (from 
Voyager UVS data, in Keel 1998). Still, the available data point to interstellar absorption 
in LBGs which is stronger and broader than in the most nearly comparable nearby systems, 
and often blueshifted compared to photospheric absorption or nebular emission lines. 

The implied abundance levels are given in Tables 4 and 5, again for the Ly a and 
luminosity groupings. The C and N abundances are roughly half solar, with no evidence 
of a trend with Lyman a strength. In this respect, the luminosity groupings may be more 
useful, since the numbers show a signiflcant increase in the C lines with luminosity. There is 
an increase in the upper hmit for N as well, with the hmit set by the fact that N II seems to 
have a floor for [0/H] below 0.48 as well as the hmited sensitivity for the fainter groupings. 

This is a milder C/0 deflcit than inferred from mid-UV emission lines by Shapley et 
al. (2003). In contrast to the C and N results, the P abundance is too low to measure (less 
than 0.1 solar) except possibly in the reddest and weakest-emission LBGs (group 1) and 
the intermediate-luminosity composite. This lack of absorption is observationally robust, 
possibly even exacerbated by the role of interstellar absorption, and suggests that these 
objects are in an early stage of nucleosynthetic history as far as this element is concerned. 
Phosphorus is a product of the reaction network in red-giant interiors, essentially from the 
combination of a- and s-processes, with some contribution from supernovae (dominated by 
type II explosions), as summarized by Trimble (1991). The high oxygen masses expected for 
SN II remnants might lead to 0/H measures that are higher than would be expected from 
scahng other elements as well, during the interval from the onset of SN II enrichment until 
the evolution of solar mass stars has time to approach the enrichment that our vicinity leads 



-9- 



us to call normal. 

Other evidence has suggested that the abundances in LBGs are low for luminous galax- 
ies, but not extremely low. In addition to the C IV wind data of Mehlert et al. (2002), 
the traditional optical bright-line technique has been applied via near-infrared spectroscopy. 
Pettini et al. (2001) showed that the gas-phase 0/H ratio must typically he above 0.1, from 

the [O III]/H/9 ratios alone. For a few selected objects with more extensive spectral coverage, 
Kobulnicky & Koo (2000) had found 0/H = 0.25-0.95 solar, fitting with the more limited 
results of Teplitz et al. (2000) from a combination of spectroscopy and fortuitously placed 
filter imaging. These results for LBGs stand in marked contrast to the low abundances found 
for damped Lyman a absorption systems, indicating a marked divergence in the history of 
these disparate environments. They also contrast with the fainter narrow Lyman a emitters 
for which colors and limited spectroscopic data suggest still lower abundances, 0/H below 
1/3 solar and conceivably as low as in I Zw 18 (Keel et al. 2002). 



4. Molecular Hydrogen Absorption 

The far-UV range offers a unique way to probe molecular gas in these young systems, 
sensitive even to cold H2 and biased only by the requirement that sufficient UV continuum 
penetrates the molecular medium to contribute to the emerging UV flux. This bias may 
account for the low limits to H2 column density derived for the M33 H II regions by Bluhm 
et al. (2003) and for I Zw 18 by Aloisi et al. (2003). In each case, analysis of FUSE data 
indicates N(H2) < 10^^ cm^^ in all cases, and < 10^^ for I Zw 18, in which the wavelength 
overlap of H2 multiplets from the target galaxy with other features in Milky Way material 
is less serious. Hoopes et al. (2001) suggest that similar column densities apply to more 
luminous starburst systems, and implicate the UV-transparency bias. This bias should be 
minimized at low abundances, where the optical and UV extinctions are smaller for a given 
gas column density (Testor & Lortet 1987; Quillen & Yukita 2001). The persistence, or 
absence, of such a bias could trace the development of the gas-to-dust ratio as elements 
needed for grain formation are progressively enriched in the interstellar medium. 

The group composite spectra of Lyman-break galaxies can be examined for possible 
signatures of H2 absorption. At their resolutions, even the strongest individual H2 features 
are blended, so I consider what combinations of H2-free template spectra (from the low- 
redshift objects and theoretical predictions) fit the LBG spectra most closely after folding 
in H2 absorption spectra for various column densities. Single column density values may 
not refiect the average across the galaxy closely, given the complex weighting of both the 
starlight and has distributions (as discussed by Bluhm et al. 2003) and the possibility that 
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actively star-forming galaxies might have multiple velocity systems which strengthen the 
overall absorption. Still, even crude values are useful given the orders-of-magnitude disparity 
between typical Milky Way sight lines and values derived for bright star-forming regions in 
nearby galaxies. The "local" templates come from the spectra of NGC 604 and I Zw 18 
after removal of galactic H2 absorption. The removal was done by fitting and subtracting 
Gaussian profiles (for weak isolated lines) or Voigt profiles (with the same parameters for 
members of each strong multiplet) at the wavelengths of all individually detected H2 lines 
from the compilation by McCandhss (2003) , except those lines confused with expected stellar- 
wind features. The ratio of spectra before and after this process gives the spectrum of 
foreground Galactic absorption sampled along these lines of sight, which we can use as 
empirical templates for H2 absorption. I also consider purely theoretical templates, using 
the absorption spectra from McCandliss (2003) for intrinsic velocity width 6= 10 and 20 km 
s~^. The H2OOLS IDL routines generate such templates for various sets of level populations, 
showing that the levels J = — 2 are well populated in the galactic foreground and would 
likely be so in other galaxies as well. Populations with J > 2, although much less populated 
than the lower levels, are generally seen to be populated well above the expected thermal 
levels based on the lower states, which is attributable to radiative excitation from the ambient 
starlight. 

Some features in the composite LBG spectra do fall near H2 bands, although the op- 
portunities for confusion with relatively weak atomic lines are many (Fig. 4). The matches 
change with assumed H2 rotational temperature, since adjacent bands shift slightly to the 
red for transitions arising at higher J. This effect amounts to typically 6 A between J = 0, 3, 
and the higher J absorption has a more complex profile (which will be broader even at the 
relatively coarse resolution of the composite spectra). This structure is one of the reasons 
that absorption from these higher-lying states is stronger for a given column density. The 
strongest constraints come from wavelengths with no observed absorption, where H2 absorp- 
tion models predict features comparably deep to the ones observed. The most prominent of 
these are near 1002, 1050, and 1078 A . On the other hand, there are features matching some 
of the strongest H2 bands in otherwise clean spectral regions near 986, 1011, 1064, and 1096 
A . Weak interstellar fines, not seen in NGC 604, do occur near these features as follows: O 
I 988, S III 1012, Si IV (wind) 1062, Fe II 1063, and Fe II 1096. The band structure leads to 
the features near 987 and 1012 becoming relatively stronger for warmer H2. This leaves two 
bands near 1050 and 1092 Aas the most sensitive regions for seeking H2 absorption. The 
highest signal-to-noise here occurs in the group 1 Ly a composite, with equivalent-width 
limits for these features of 0.27 and 0.20 Arespectively. 

For the group 1 composite spectrum, comprising the LBGs with net absorption at 
Lyman a, regions of the spectrum not affected by strong atomic lines can be matched with 
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a synthetic spectrum for N(H2) — 1.7 x 10^^ cm~^ if the molecules are mostly in J = 3; this 
fit does not violate any spectral regions with high local S/N. This J comes from the mean 
wavelengths of spectral dips near 1098 and 1114 A , as well as the strength of the 987 and 
1012-A features. These data rule out comparable contributions from lower states at more 
stringent equivalent-width limits for decreasing J. However, the limits on N(H2) in the lower 
states are comparable, since the smoothed absorption per unit column density is weaker for 
these transitions. 

For the group- 1 composite, one could view this as evidence for a substantial amount of 
warm molecular absorption, or as an upper limit including lower states of ~ 7 x 10^^ cm^. 
This is a very modest column density by the standards of Milky Way sight lines; the FUSE 
spectra of I Zw 18 and the M33 H II star- forming regions both show foreground absorption 
at N(H2) = 5 X 10^^ cm ~^ at galactic latitudes b = 44° and -31° respectively. In both 
of these objects, low hmits to internal H2 absorption have been set from FUSE data (Aloisi 
et al. 2003, Bluhm et al. 2003, Keel et al. 2004). The cleanest mcasTu^cment is < lO^'' cm'^ 
for I Zw 18, at a large enough redshift to reduce contamination by Milky Way absorption. 
Limits for the M33 regions are < 10^'' cm~^ in all cases. All these values must be affected by 
the UV bias; the mixing of molecular gas with dust means that stars behind large amounts 
of H2 will not contribute to the overall UV flux. 

The other groups of LBGs have slightly lower limits on H2 absorption than the possible 
detection in the group 1 objects. This result indicates that the ratio of UV extinction to 
molecular-gas column density (that is, the dust-to-gas ratio) was not much smaller in the 
Lyman-break galaxies than in local star-forming systems, so that stars behind substantial 
amounts of H2 are likewise obscured and vanish from the UV spectrum as well. This in turn 
suggests that any large amounts of H2 in Lyman-break galaxies are mixed with dust (and 
by implication, the whole range of metals) rather than being primordial and unenriched, 
because it is just the latter case which would show up most effectively via strong absorption 
in the Lyman and Werner bands. 

The Lyman a forest will limit our ability to derive more precise H2 column densities, 
since in any individual object at high redshift, there will be confusion between its absorption 
and random H I features. We may be limited to such rough "typical" class properties. 

5. Conclusions 

Far-UV spectra of nearby star-forming galaxies are compared to composite spectra of 
Lyman-break galaxies at 2; ~ 3, with special attention to abundance-sensitive features. Fea- 
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tures of C and N suggest moderately subsolar abundances, at the level associated with 0/H 
roughly half solar in present-day galaxies. There is some evidence that the C abundances 
correlate with mid-UV luminosity, and N may do so subject to detection limitations that 
introduce upper limits. Phosphorus is underabundant by an order of magnitude in most 
subsamples of LBGs, perhaps reflecting a dominant origin in red giants rather than SN II. 
This fits with additional evidence from emission-line measurements, that the ionized gas is 
of low but not very low abundances, and suggests in itself that the star-forming histories of 
LBGs have been brief but eventful. 

The composite spectra also allow limits on the column density of H2 against the UV 
starlight in typical LBGs, which are comparable to what we see looking out of the Milky Way 
at modest latitudes I ^ 30°. While some bias toward unextinguished stars must exist, this 
suggests that there are not huge reservoirs of cold and unenriched gas surrounding LBGs. 
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LBG spectra in digital form, and for extensive and challenging discussions. This work 
was supported in part by NASA FUSE grant NAG5-8959, with pubhcation support kindly 
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service rendered by Steven McCandliss in making his H2OOLS software and data compilation 
publicly accessible. Alessandra Aloisi provided the summed I Zw 18 spectra, and B.-G. 
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Fig. 1. — A field guide to far-ultraviolet galaxy spectra. The FUSE spectra of NGC 604 
(thick line) and I Zw 18 (thin fine) are shown, with Galactic H2 absorption fines fitted 
and removed, and smoothed by 0.3 A for graphic simplicity, with various line wavelengths 
illustrated. Wavelengths of the strongest H2 transitions are shown by the lowest set of marks; 
the line list is the set that was strong enough for Milky Way absorption to be detected against 
NGC 604. The strongest atomic absorption features from the interstellar medium are marked 
in the middle row. Lines with significant stellar-wind contributions are marked at the top 
with schematic blueward wings on each symbol, and ions identified. Lyman {3 lies at the 
break between panels. The relative fiux is in F^, and represents the observed fiux in erg cm~2 
s-^ A scaled upward by 10^^ for NGC 604 and 2 x lO^'' for I Zw 18. Differences between the 
two spectra include stronger stellar- wind features, in both emission and absorption portions 
of the P Cygni profiles, in NGC 604, stronger absorption in P V from NGC 604 for which the 
dominant contributor is stellar photospheres, and a few apparently blueshifted interstellar 
features in I Zw 18 which are foreground Milky Way gas. The wavelength scale is the emitted 
frame for each object. 
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Fig. 2. — Comparison of the composite LBG spectra with the FUSE data on NGC 604 and 
I Zw 18, after foreground atomic and H2 absorption hnes have been removed. The LBG 
groups are numbered as in Shapley et aL (2003) and Table 2, with group 1 having net 
Lyman a absorption, and 2-4 with stronger net emission in that order. All the spectra have 
been boxcar-smoothed to a common resolution 4.1 A for display, a value set by convolving 
the LBG composites with their typical resolution to enhance the visibility of weak lines. 
The flux zero point is for the lowest spectrum; successively higher ones are offset by integer 
values units. Important interstellar and photospheric features are marked at the bottom; 
wind lines are shown across the top. The symbols accompanying the wind designations are 
intended to suggest a blueshifted absorption trough, but some of these lines have distinct 
redshifted emission components as well (often blended with components of additional lines, 
as in the O VI doublet). 
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Fig. 3. — Comparison of the composite LBG spectra, now grouped in three bins of mid-UV 
luminosity, with the FUSE data on NGC 604 and I Zw 18, displayed as in Fig. 2. Several 
strong lines show systematic trends with UV luminosity. For example, the O VI doublet, 
N II A1083, and N III A991 all strengthen for brighter galaxies,. Brighter galaxies are also 
redder in continuum slope shortward of about 1000 A. 
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Fig. 4. — Composite LBG spectra, with the continuum shape normahzed to indicate absorp- 
tion depths, compared to theoretical absorption spectra, as calculated with S. MacCandliss' 
H2OOLS routines. Each J value is shown separately, for a column density 1.7 x 10^^ cm~^. 
This value is the maximum that any J population could have against the group 1 spec- 
trum without violating the data, although having populations weighting strongly to J = 3 
is physically implausible. Potentially confusing ionic absorption wavelengths are labelled at 
the bottom. 
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Table 1: Prominent Far-UV Features in Galaxy Spectra 



Note. — Wavelengths are as emitted. The dominant source of each hne is given as photospheric, winds, 
or ISM, following the attributions in Leitherer et al. (2002), Walborn et al. (2002), and Danforth et al. 
(2002). A superscript 3 denotes that the listed wavelength is for the central member of a triplet, and an 
asterisk indicates absorption from an excited state. 
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Table 1: Far-UV Absorption Features in Emission- Grouped Composite LBG Spectra 



Note. — Equivalent widths in A in the emitted frame; positive values represent absorption. 



Group: 
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Table 2: Far-UV Absorption Features in Luminosity-Grouped Composite LBG Spectra 



Note. — Equivalent widths in A in the emitted frame; positive values represent absorption. 
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Table 3: Derived Effective [0/H] for Emission-Grouped Composite LBG Spectra 



Note. — Abundances are in solar units, with a local calibration referred to nebular [O /H] measurements. 
Multiple ions are combined by averaging in the bottom lines. 
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Table 4: Derived Effective [0/H] for Luminosity- Grouped Composite LBG Spectra 



Note. — Abundances are in solar units, with a local calibration referred to nebular [0/H] measurements. 
Multiple ions are combined by averaging in the bottom lines. 



